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The a-proton Coulomb interaction is included in the description of a-deuteron scattering using 
the screening and renormalization approach in the framework of momentum-space three-particle 
equations. The technical reliability of the method is demonstrated. Large Coulomb-force effects are 
found. 
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I. INTRODUCTION 

The application of exact Faddeev three-body theory 
to the understanding of nuclear reactions that are dom- 
inated, at low energies, by three-body degrees of free- 
dom has been shadowed in the past by the difficulty in 
dealing with the long-range Coulomb interaction between 
charged particles. Up to now this difficulty has been 
overcome by the use of approximate theoretical meth- 
ods, namely, the continuum-discretized coupled-channel 
(CDCC) method El, 

which was proposed initially to deal 
with deuteron scattering on a heavier nuclei, but that is 
now widely used to analyze data resulting from the reac- 
tions involving halo nuclei. Given the progress achieved 
recently in the description of proton-deuteron (pd) elas- 
tic scattering and breakup using exact three-body equa- 
tions 0, EJ 0] , we are now able to include the Coulomb 
interaction in three-body nuclear reactions involving two 
charged particles. Before we apply the present method to 
the description of direct nuclear reactions that are dom- 
inated by three-body degrees of freedom, we show here 
the results for alpha-deuteron (ad) scattering which at 
low energies is one of simplest effective three-body nu- 
clear reactions. 

The ad scattering has been studied extensively in the 
past, both experimentally and theoretically H, 0, H, 
SEE E1II1- Although' ad is a six-nucleon system, at 
low energies, to a good approximation, the a particle 
may be considered a spin zero structureless boson, and 
thereby the theoretical description of ad scattering may 
be reduced to a three-body problem made up of one a 
and two nucleons (N). Then, the most serious difficulty 
is the treatment of the long-range Coulomb interaction, 
which in the previous calculations was either completely 
neglected or taken into account only approximately. Us- 
ing the method developed in Refs. 0, y, H[ we are now 
able to include the Coulomb interaction also for ad scat- 
tering, quantitatively evaluate its importance, and with 
greater confidence ascertain the quality of aN force mod- 
els. 

Section |TT] shortly recalls the technical apparatus un- 
derlying the calculations. Section [TIT] presents character- 



istic results. Section ITVl gives our conclusions. 



II. TREATMENT OF THE COULOMB 
INTERACTION USING THE SCREENING AND 
RENORMALIZATION APPROACH 

Our treatment of the Coulomb interaction is based 
on the idea of screening and renormalization proposed 
in Ref. 1X3 1 fo r two-particle scattering and extended in 
Refs. @,M1I3>EB| to three-particle scattering where only 
two particles are charged. The Coulomb potential is 
screened, standard scattering theory is applicable, and 
the renormalization procedure is applied to recover the 
unscreened limit. The success of the method depends on 
the choice of the screened Coulomb potential 



w R (r) = w(r) 



-(r/Rf 



(1) 



We prefer to work with a sharper screening than the 
Yukawa screening (n = 1) of Ref. EH . We want to ensure 
that the screened Coulomb potential wr{t) approximates 
well the true Coulomb one w(r) for distances r smaller 
than the screening radius R and simultaneously vanishes 
rapidly for r > R, providing a comparatively fast con- 
vergence of the partial-wave expansion. In contrast, the 
sharp cutoff [n — > oo) yields an unpleasant oscillatory 
behavior in the momentum-space representation, leading 
to convergence problems. In Refs. 0, Q we found the 
values 3 < n < 6 to provide a sufficiently smooth, but 
at the same time a sufficiently rapid screening around 
r = R; n = 4 is our choice also in the present paper. 

We solve Alt-Grassberger-Sandhas (AGS) three- 
particle scattering equations Ell i n momentum space 

U^\z) = S 0a G^(Z) + J2h^ R} (Z)Go(Z)U^(Z) : 

(2a) 

U£\Z) = Gv\Z) + J2T^(Z)G (Z)U^(Z), (2b) 



with 5/} a = 1 
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5/3q, Gq(Z) being the free resolvent, 
To (Z) the two-particle transition matrix derived from 
nuclear plus screened Coulomb potentials, and uj^\Z) 

( R) 

and Uqq, (Z) the three-particle transition operators for 
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elastic/rearrangement and breakup scattering; their de- 
pendence on the screening radius R is notationally in- 
dicated. On-shell matrix elements of the operators ([2]) 
between two- and- three-body channel states \4>a{(\i)v ai ) 
and |0o(p/q/)^O/} with discrete quantum numbers v aj , 
Jacobi momenta pj and cjj, energy Ei, and Z = Ei + iO, 
do not have a R — > oo limit. However, as demonstrated in 
Refs. [2j, |3j, llJ] , the three-particle amplitudes can be de- 
composed into long-range and Coulomb-distorted short- 
range parts, where the quantities diverging in that limit 
are of two-body nature, i.e., the on-shell transition ma- 
trix T£ii(Z) derived from the screened Coulomb poten- 
tial between spectator and the center of mass (cm.) of 
the bound pair, the corresponding wave function, and the 
screened Coulomb wave function for the relative motion 
of two charged particles in the final state. Those quanti- 
ties, renormalized according to Ref. [13], in the R — > 00 
limit converge to the two-body Coulomb scattering am- 
plitude (<£a(q/)f Q/ \T"c\<l> a (<±i)v ai ) (in general, as a dis- 
tribution) and to the corresponding Coulomb wave func- 
tions, respectively, thereby yielding the three-particle 
scattering amplitudes in the proper Coulomb limit 



(<A/3( t l/) z/ / 8/ \Up a \4>a{<li)Vai) 

= 8j3 a (<f> a (c[f)v af \ T£c~ 1 4> a (q< ) v ai ) 

+ £?J Z PR (?/) to (<i/K I Ha ( E * + *o) 

(3a) 

{MPfClf)vOf\Uo a \<l> a (qi)v ai ) 

_i 

= D lim { z r* (Pf)iMPf<lf) v o f \ (3b) 

ft — »00 v ' 

x U^ R) (Ei +{0)^(^)^)2 J( Qi )}. 

The renormalization factors Z a n(qj) and ZR(p f) are di- 
verging phase factors given in Refs. 0, [H, [llj 13 ■ The 
R ! — ► 00 limit in Eqs. ([3]) has to be calculated numer- 
ically, but due to the short-range nature of the corre- 
sponding operators it is reached with sufficient accu- 
racy at rather modest R if the form of the screened 
Coulomb potential has been chosen successfully as dis- 
cussed above. More details on the practical implementa- 
tion of the screening and renormalization approach are 
given in Refs. [1, H, Compared to the pd calculations 
an additional difficulty is the presence of the aN P-wave 
transition matrix T„(Z) poles in the complex energy 
plane close to the real axis. This requires special treat- 
ment of those partial waves in the numerical solution of 
Eqs. @, i.e., the pole factor olT„ R \z) is separated when 
interpolating Ta R \z)Go(Z)Ua R \z) and the subtraction 
technique is used for integration. 



III. RESULTS 

There exist in the literature several parametrizations 
of the aN potential. Most of them support deeply bound 
aN state in S-wave that is not observed experimentally, 
since it is forbidden by the Pauli principle. To account 
for the Pauli principle that forbidden state \b) (or some 
state close to it) has to be projected out or moved to 
a large positive energy T, replacing the potential v by 
v' = v + |fr)r(6|. As demonstrated in Ref. [HI], the lat- 
ter method in the r — > 00 limit is equivalent to the first 
one. Alternatively, one may use repulsive aN potential in 
the S'-wave. In order to estimate the model dependence, 
in the calculations of this paper we use three different 
parametrizations for the nuclear part of the aN poten- 
tial, that in the following will be called aN-l, aA^-II, and 
aN-111. The potential aN-l is taken from Ref. [l8]. It is 
of Woods-Saxon form with central and spin-orbit parts 
and supports a bound-state in S'-wave that is moved to a 
large positive energy T; we found that for T > 1000 MeV 
the results are practically independent of T. The poten- 
tial ctA^-II is taken from Ref. [18[ as well and differs from 
aN-l only in S-wave where it is repulsive and local, but, 
nevertheless, nearly phase equivalent to aN-l. The third 
potential aA^-III, taken from Ref. [l(|, is represented as 
a sum of Gaussians with the strong repulsive rank-1 sep- 
arable term in S-wave that moves the harmonic oscil- 
lator ground state which is close to the Pauli-forbidden 
state to T = 1000 MeV. The potentials fit the low en- 
ergy experimental aN phase shifts quite well, although 
aiV-III is more attractive in P-waves than the other two. 
All potentials are charge symmetric and act in partial 
waves with orbital angular momentum L < 2. For the 
neutron-proton (rip) interaction we use the charge depen- 
dent (CD) Bonn potential [lj| and include partial waves 
with total angular momentum I < 2 plus 3 D^. It was 
checked that those quite low partial waves are sufficient 
for convergence. In contrast, the screened Coulomb po- 
tential is of longer range and therefore, depending on the 
screening radius R, requires partial waves up to L < 11. 
Those high partial waves are included exactly as de- 
scribed in Ref. [1| . We note that the ap Coulomb poten- 
tial at short distances is taken as the one of a uniformly 
charged sphere, i.e., w(r) = 2a e [3 — (r/Rj) 2 ]/(2Rj) for 
r < Ri — 1.5 fm, a e being the fine structure constant. 

The internal criterion for the reliability of the screen- 
ing and renormalization approach is the convergence of 
the observables with the screening radius R employed 
to calculate the Coulomb-distorted short-range part of 
the amplitudes in Eqs. ([3]); that criterion was found to 
be absolutely reliable for pd scattering Q • Figures Q] 
and [2] show several examples for ad elastic scattering and 
breakup. The breakup kinematical final-state configura- 
tions are characterized in a standard way by the polar an- 
gles of the a-particle and the proton and by the azimuthal 
angle between them, (9 a , 9 P , tp p — tp a ). The definition of 
the arclength S alongthe kinematical curve is the stan- 
dard one as in Refs. [1, [l(J. Even when the Coulomb 
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FIG. 1: Convergence of the ad elastic scattering observables 
with screening radius R. The differential cross section and 
deuteron analyzing powers iTn and T20 at 4.81 MeV deuteron 
lab energy are shown as functions of the cm. scattering angle. 
Results obtained with screening radius R — 5 fm (dotted 
curves), 10 fm (dash-dotted curves), and 15 fm (solid curves) 
are compared. Results without Coulomb (dashed curves) are 
given as reference for the size of the Coulomb effect. The 
calculations were performed with the aiV-I potential. 



effect is large the convergence is impressively fast, e.g., 
the screening radius R — 10 fm is practically sufficient 
for the elastic ad scattering at 4.81 MeV deuteron lab 
energy in Fig. [TJ while the observables of ad breakup at 
15 MeV a lab energy in Fig. [2] require at least R = 15 fm 
for convergence. As discussed in Refs. @, 0|, the con- 
vergence rate is energy dependent, i.e., larger screening 
radii are needed for elastic scattering observables at very 
low energies and for the breakup differential cross sec- 
tion in kinematical situations characterized by very low 
relative energy of the two charged particles, i.e., close to 
the ap final-state interaction (FSI) regime. Nevertheless, 
the observed convergence strongly suggests the reliabil- 
ity of the present Coulomb treatment using screening and 
renormalization approach. 

As already shown in Figs. [T] and [2] the Coulomb effect 
on the observables of low energy ad scattering may be 
very strong, indicating that the inclusion of the Coulomb 
interaction is necessary for a stringent comparison of the 
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FIG. 2: Convergence of the ad breakup observables with 
screening radius R. The differential cross section and 
deuteron analyzing power A yy in selected kinematical con- 
figurations at 15 MeV a lab energy are shown as functions 
of the arclength S along the kinematical curve. Results ob- 
tained with screening radius R — 10 fm (dotted curves), 15 fm 
(dash-dotted curves), and 20 fm (solid curves) are compared. 
Results without Coulomb (dashed curves) are given as refer- 
ence for the size of the Coulomb effect. The calculations were 
performed with the aN-1 potential. 



theoretical results and experimental data. Obviously, we 
have many more predictions than it is possible to show. 
Therefore we make a selection of the most interesting 
predictions which illustrate the message we believe the 
results tell us. The readers are welcome to obtain the 
results for their favorite data from us. 

Figure [3] presents our results for the differential cross 
section and all deuteron analyzing powers of elastic ad 
scattering at 4.81 MeV and 17 MeV deuteron lab en- 
ergy. As expected, the Coulomb effect is large at lower 
energy for all observables in the whole angular regime, 
while with increasing energy it becomes smaller. Fur- 
thermore, the results depend strongly on the choice of 
the aN interaction. The potential aN-1 describes most 
of the experimental Ed — 4.81 MeV data quite satisfac- 
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FIG. 3: (Color online) The differential cross section and deuteron analyzing powers at 4.81 MeV and 17 MeV deuteron lab 
energy. Results including the Coulomb interactions obtained with potentials aN-1 (solid curves), aiV-II (dotted curves), and 
QfTV-III (dash-dotted curves) are compared. aN-I results without Coulomb (dashed curves) are given as reference for the size 
of the Coulomb effect. The experimental 4.81-MeV data are from Ref. [3] for the differential cross section and from Ref. 0] for 
the spin observables, and the 17-MeV data are from Ref. [f|. 



torily, whereas the potential aiV-II that avoids the Pauli 
forbidden state in S-wave by the local repulsion clearly 
fails in accounting for the scattering data, especially at 
Ed = 17 MeV, although the 6 Li bound state properties 



predicted by both potentials are very similar. The poten- 
tial aN-lll also fails in accounting for the experimental 
data at Ed — 4.81 MeV, but the difference between pre- 
dictions using aN-I and aN-III is found to be mostly due 
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to the differences in P-waves. At Ed = 17 MeV, which 
is just 3 MeV below 3 H + 3 He threshold, there is a quali- 
tative agreement between experimental data and predic- 
tions with aN-1 and aJV-III, except for deuteron tensor 
analyzing powers T20 and T21. Obviously, the treatment 
of the a-particle as a structureless boson becomes less 
reliable with increasing energy. 

Figure H] presents our results for the differential cross 
section of ad breakup at 15 MeV a lab energy in var- 
ious kinematical configurations. The most important 
Coulomb effect is the shift of the ap P-wave resonance 
position that leads to the corresponding changes in the 
structure of the observables. Furthermore, the Coulomb 
interaction breaks an — ap charge symmetry and thereby 
allows the coupling to the np isospin triplet waves, in par- 
ticular So- The predictions without Coulomb fail com- 
pletely in accounting for the experimental data, while 
inclusion of the Coulomb moves the peaks of the differ- 
ential cross section to the right positions, although the 
height of those peaks is not always reproduced. There 
is a strong dependence on the employed aN potential, 
where, like in the elastic scattering, the potential aiV-II 
provides the worst description of the data, and the dif- 
ference between the predictions of aN-1 and aTV-III is 
dominated by the P-waves. The sensitivity of the results 
to the choice of a realistic np potential is also checked 
and found to be insignificant. When the cm. energy in- 
creases one finds a larger part of the phase space where 
the relative ap energy may be quite different from the 
one corresponding to the P-wave resonance and there- 
fore also the Coulomb effect is less significant as shown 
in Fig. [5l 

IV. SUMMARY 

In this paper we show that the screening and renor- 
malization approach for the inclusion of the Coulomb in- 



teraction in the description of three-particle scattering 
using momentum-space integral equations, developed in 
Refs. @,H,IU f° r pd reactions, can be extended reliably to 
ad scattering. This is an important step towards appli- 
cation of exact scattering equations for the description of 
nuclear reactions within three-body models, where up to 
now only approximate treatments like the CDCC method 
m have been used. 

The Coulomb effect on the observables of ad scattering 
is studied and is found to be large in elastic scattering 
at very low energies and in breakup, where the shift of 
ap P-wave resonance position leads to the corresponding 
shifts of the differential cross section peaks. 

Another important consequence of this work is that it 
allows to ascertain with greater confidence the quality 
of the aN force models one uses to describe ad observ- 
ables as well as structure and reactions of 6 Li and 6 Hc. 
Although at present there are too large uncertainties in 
the parametrization of P- and P-wave interactions that 
need to be improved, this work clearly indicates the su- 
periority of the attractive S'-wave potentials supporting 
a Pauli-forbidden state (that is projected out or moved 
to a large positive energy) over the repulsive .S'-wave po- 
tentials. 
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FIG. 4: (Color online) The differential cross section of the ad breakup at 15 MeV a lab energy in selected kinematical 
configurations. Curves as in Fig. [3] The experimental data are from Ref. [8j. 
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FIG. 5: The differential cross section and deuteron analyzing power A xx of the ad breakup at 12 MeV deuteron lab energy 
in the (25°, 45°, 180°) configuration. Results obtained using the potential aN-1 with (without) Coulomb are shown as solid 
(dashed) curves. The experimental data are from Ref. 0. 
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